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ABSTRACT Singlet oxygen generated by optical pumping in 
liquid oxygen/air medium has recently been reported as a po¬ 
tential gain medium for high power lasers at 1580 nm, by 
USA’s defense agency DARPA. However, the details with ref¬ 
erence to the underlying physics of this laser and the potential 
pumping techniques for achieving lasing are still unclear. The 
present paper investigates numerically the population kinetics 
of both upper and lower lasing levels viz. O2 ( 1 A g ) i;= o and 
02( 3 U g ) v =i states in liquid oxygen and discusses the criticality 
involved in obtaining lasing with this medium. Isotopic liquid 
oxygen and liquid air medium where improved conditions for 
lasing are anticipated, as compared to that in natural liquid oxy¬ 
gen, have also been considered in the present study. The studies 
have been carried out for optical pumping by both continuous 
and pulsed mode lasers at 1064 nm and 634 nm wavelengths. 
The temporal variation of small signal gain in each case has 
been investigated along with limitations of the same from point 
of view of lasing. The available and extractable power from all 
three medium for an optimum pumping case has also been dis¬ 
cussed. The studies reveal that liquid air has high potential for 
very high power lasers as compared to both natural and isotopic 
liquid oxygen. 
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1 Introduction 

The search for an ideal lasing species with en¬ 
ergy levels having appropriate meta-stable lifetime has been 
a quest for researchers from the day of laser invention. High 
power continuous wave lasers (multi-kilowatt) have wide ap¬ 
plications in both civil as well as defense areas. Even though 
megawatt class continuous wave lasers: hydrogen/deuterium 
fluoride lasers (HF/DF), chemical oxygen iodine lasers 
(COIF), CO2 gas dynamic laser (GDF), etc. are available at 
present, their applications to a large extent have been limited 
by issues pertaining to compactness. Further, handling of haz¬ 
ardous chemicals, operational easiness and readiness time 
are the other practical issues related to the above high power 
lasers. These limit their applications and make it imperative to 
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search for new lasing mediums. Therefore, the search for an 
extremely low specific weight (weight per unit power output) 
and eco-friendly laser still remains a front line area of interest 
among laser researchers. The first electronically excited level 
of oxygen molecules, the singlet delta oxygen O2 ( ! Ag), is 
one such species which fulfills the desired properties of a near 
ideal lasing medium. Various methods such as optical, elec¬ 
trical and chemical have been utilized for the generation of 
these singlet delta oxygen molecules. Among these, the most 
efficient method is chemical one and generated singlet delta 
oxygen has been efficiently used as the pumping medium in 
COIF [1]. However, all these singlet oxygen generation tech¬ 
niques have failed to yield a stimulated emission condition, 
in either of its strong emission lines at 1268 nm or 1580 nm. 
Until very recently, the research in singlet oxygen generation 
had primarily advanced from the viewpoint of its application 
as a pumping source in COIF rather than as a lasing medium. 

The exploration of liquid oxygen as a lasing medium 
started in mid-1970s [2]. But success in terms of lasing was 
reported only recently, in 2005, by USA’s DARPA organi¬ 
zation [3]. The physics and techniques of this laser have 
not been published. Optical excitation of the gas phase oxy¬ 
gen molecules is extremely inefficient due to its very low 
absorption coefficient. In the liquid phase (97 K) the oxy¬ 
gen molecules exist in ‘pair’ form, which increases their ab¬ 
sorption coefficient many folds for its various characteristic 
lines [4]. This makes optical pumping in liquid oxygen effi¬ 
cient and hence opens up the possibility of very high density 
production of meta-stable oxygen molecules. During such op¬ 
tical pumping, successive electronic to vibrational transfer 
(E-V transfer) and vibrational to vibrational (V-V) relaxation 
of the higher levels is found to occur. The C>2( 3 U g ) v =i is an¬ 
other meta-stable level which has a long lifetime of ~ 2.5 ms 
in case of pure natural liquid oxygen. This leads to accu¬ 
mulation at the longer living first vibrational level of the 
ground molecules, C>2( 3 ZJ g ) v =i- Hence, under steady state 
condition, liquid oxygen medium gives way for a three-level 
laser configuration with emission between 02( 1 A g ) i;= o and 
02( 3 Z} g ) v =\ levels at 1580 nm. 

In this paper, detailed investigation on the population ki¬ 
netics of the upper [02( 1 A^) u= o] and lower [02( 3 X , g) u =i] 
lasing levels have been presented. The procedure employed 
considers various decay processes along with setting up the 
corresponding rate equations utilizing the observed rate con- 
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stants by various experimental works. Continuous as well as 
pulsed pumping options employing sources having 1064 nm 
and 634 nm wavelengths with various power densities are 
considered. Further, the analysis has been conducted for three 
different oxygen rich cryo-mediums considering their corres¬ 
ponding kinetics to arrive at the best possible option. The 
temporal evolutions of the small signal gain in each case and 
power extraction analysis for the best possible cases are also 
discussed. 

2 Laser chemical kinetics and physics 

The absorption and emission lines of liquid oxy¬ 
gen have been explained theoretically and experimentally by 
many researchers [2,4,5]. Detailed study on the absorption 
lines of oxygen in the liquid phase and gas at very high 
densities is reported by Zhdanov et al. [4]. The above stud¬ 
ies clearly indicate that the absorption cross section is im¬ 
proved many folds for almost all lines in liquid oxygen as 
compared to that in gaseous oxygen, which confirms the ex¬ 
istence of pair molecules in liquid phase. Considering the 
existence of pair molecules in liquid oxygen, major energy 
level scheme for co-operative emissions and absorptions in li¬ 
quid oxygen is shown in Fig. 1. The major absorption lines of 
liquid oxygen are 470 nm, 581 nm, 634 nm, 764 nm, 1064 nm 
and 1268 nm with relative absorption coefficients being 2.4, 
7.1, 5.1, 1.0, 5.4 and 6.8 respectively [4]. On the other hand 
the major emission lines are the following: 581 nm, 634 nm, 
703 nm, 1268 nm and 1580 nm, having relative emission in¬ 
tensities [2,6] as 0.03, 1.0, 1.0, 3.0 x 10 4 and 1.0 x 10 4 re¬ 
spectively. This indicates that the emissions at 1268 nm and 


1580 nm form the maximum probable emission lines in the li¬ 
quid oxygen medium from viewpoint of lasing. The strongest 
emission line, 1268 nm, is also a strong absorption line in 
liquid oxygen limiting its chance of being a lasing line. There¬ 
fore, 1580 nm emission is expected to be the prominent lasing 
line as it has almost zero absorption at this wavelength. 

In liquid oxygen, the lifetime of the higher electronic state 
(02 1 27), the vibrational levels of the first electronic level, 
(O 2 l A) v = 1 - 3 , and the higher vibrational levels of the ground 
state molecules, (O 2 3 £) v> u are very low and hence quench 
and relax very rapidly. This is well supported by various ex¬ 
perimental studies [7,8], where no spontaneous anti-Stokes 
Raman scattering lines from any vibrational levels (O 2 3 £) v >i 
were observed. The only line observed was for O 2 3 £ v =i level. 
Therefore, under steady state condition, simply a three level 
scheme in liquid oxygen can be considered as shown in Fig. 2. 

Correspondingly, the stimulated emission process can be 
expressed as: 

0 2 (A)„ =0 + l/iu(1580nm) 0 2 ( 3 ^) u=1 +2 hv. (1) 

Considering the strong absorption lines in liquid oxygen, 
sources with wavelengths 581 nm, 634 nm, 1060 nm and 
1268 nm sources can be used for the excitation to the first ex¬ 
cited electronic level of oxygen. Hence, the pumping to the 
upper lasing level can be carried out with any of these follow¬ 
ing photon excitation [5,9] processes: 

20 2 ( 3 X , g ) l , =0 + /!u(581 nm) -» 0 2 ( 1 ^) u=0 

+ 0 2 ( l A g ) v=1 , (2) 
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FIGURE 1 Energy level scheme in liquid oxygen 
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FIGURE 2 A possible three level lasing scheme in liquid oxygen 


2O2 ( 3 X' g ) i;= o + hv (1064 nm) - 

02 (A)„ =0 



+ O2 ( 3 X g )„ = i , 

(3) 

2O2 ( 3 ^ g ) v=0 + hv( 1064 nm) - 

->o 2 (U,)„ =1 



+ O2 ( 3 r g )„ =0 , 

(4) 

20 2 ( 3 rg) u=0 + /*u(634nm) - 

2O2 ( 1 A g ) v=0 , 

(5) 

O2 ( 3 £ g ) v=0 + hv( 1268 nm) - 

O2 ( 1 A g ) v=0 . 

(6) 


In processes (2) and (4), even though one ( l A g ) v= \ level is 
generated for each photon, it will rapidly quench to the needed 
( 1 A g ) i;= o level through the process, 

O 2 ( l A g ) v=x +O 2 ( 3 ^) u=0 O 2 ( 1 A g ) v=0 + O2 ( 3 U g ) v= i , 
k = 3.4 x 10 -14 (7 o /7) - 2 cm 3 mol _1 s _1 . (7) 


The experimentally verified refractive index of the liquid oxy¬ 
gen for the visible region (4416A to 6705A) varies from 1.225 
to 1.221, and for infinite wavelength it is estimated to be 
1.22 [ 11 ]. The collision-induced line widths for the absorption 
and emission bands in oxygen have been discussed in detail 
in [12]. The line width is nearly 200 cm -1 for most pump¬ 
ing options and by considering a radiative lifetime of 4 s, the 
stimulated emission cross section for the laser emission at 
1580 nm (6329 cm -1 ) is calculated to be 2.7 x 10 -23 cm 2 . 

Once the population inversion condition is achieved, the 
available (stored) laser power within the medium will be pro¬ 
portional to the small signal gain (y) and the saturation inten¬ 
sity (/ sat ) of the medium and is expressed as 

Pav = haiYLA. (11) 


The small signal gain and saturation intensities can be esti¬ 
mated using the relations 


y(v) = otui(u) 




hv 

ax l 1 + i) 



( 12 ) 

(13) 


However, it will be disadvantageous for lasing, since it 
simultaneously populates the lower lasing level, i.e., one 
02 ( 3 £’)u=i molecule will also be generated along with the 
desired 02( 1 A g ) u= o. Therefore, overall optical pumping in li¬ 
quid oxygen is a complex process, and hence the numerical 
studies are essential to determine the population in the cor¬ 
responding energy levels considering all possible chemical 
kinetics. 

The stimulated emission is possible only when the pop¬ 
ulation inversion condition is achieved and its probability is 
determined by the stimulated emission cross section (cr u i), 
which is estimated using the relation [10], 


CTul(u) = 




Here A u i is the Einstein coefficient for the spontaneous 
emission characterized by the upper level radiative lifetime 
(r r ), n is the refractive index of the medium at the emis¬ 
sion wavelength A and g(v) is the normalized spectral line 
shape function at the emission frequency, v. Considering that 
the medium operates at cryo-temperatures, the thermal and 
Doppler broadening will not be significant. Hence it is logical 
to use the line shape of the pump beam for present estimation 
of cr u \. Therefore, the normalized Gaussian line shape function 
is considered for this estimation [10], 


g(v) 


2 / In 2 A 
Av \ 7T ) 



where the line width A v can be written as (c/Aq) A A and hence 
(9) can be modified as 


CTul(u) = 


A 3 A 1 
4cTtn 2 AA i r 


In 2 


Tt 


( 10 ) 


where Na , A^vib are the population in the upper and lower las¬ 
ing level i.e., [0 2 ( ^g)^] and [0 2 ( 3 2 ^)i;=i], respectively, 
r u is the lifetime of the upper level, g u and g\ are their cor¬ 
responding degeneracies. For most lasing mediums, the ratio 
of the degeneracies varies from 0.5 to 1, and here we consid¬ 
ered the worst scenario, which gives a slight underestimation. 
Therefore, the small signal gain and saturation intensity can 
be rewritten as 


y(v) = cr u i(u) 




~ (tf) \_N a - A^vib] , 



hv 


err u 



hv 
2 <7T U 


(14) 

(15) 


Based on this, the gain throughout this study is simply calcu¬ 
lated from the product of the difference in population densities 
between the lasing levels (dAO with the stimulated emis¬ 
sion cross-section 2.7 x 10 -23 . The saturation intensity for the 
medium is estimated to be 53 MW/cm 2 , which is almost 20 
times less than that of its breakdown intensity. The very low 
stimulated emission cross-section urges for very high pump¬ 
ing requirements, at the same time the high saturation inten¬ 
sity also yields the possibility of efficient power extraction 
from the medium. 


3 Numerical modeling-methodology 

The detailed numerical studies have been carried 
out to determine the temporal behavior of population in upper 
and lower lasing levels by considering all major source and 
decay mechanisms corresponding to each level. The govern¬ 
ing rate equation for the population in upper and lower lasing 
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levels are expressed as follows: 


dNu 

d t 
d Ni 

d t 


r-k 1 N u -k 2 N 2 u -yN u , 

(16) 

5 

r + —k\N u + yN u — k 3 N \, 

(17) 


where N u and N\ are the population of the upper 02( 1 A g ) v= o 
molecules and lower 0 2 ( 3 £ g ) v= i molecules, respectively, 
and k is the rate constant for the various probable decay or 
quenching processes in liquid oxygen: 


O2 ( 1 Z\ g ) y=0 + O2 ( 3 X'g) y= o —> 2O2 ( 3 X' g ) y= o ^ 

k\ = 1.05 x 10 -18 cm 3 mol -1 s -1 , (18) 

O2 ( 1 A g ) v=0 + O2 ( x A g ) v= Q —> 2O2 ( 3 Z g ) v=0 + 635 nm, 

_ 1 

k 2 = 1.7 x 10 -18 cm 3 mol _1 s _1 , (19) 

0 2 ( l A g ) v= 0 -> 0 2 C^ g ) v =\ + 1580nm, 

y = \ s_1 , ( 2 °) 

0 2 ( 1 ^) y=0 + 4O 2 ( 3 i: g )v=0 -> 50 2 ( 3 U g ) v= i 
+ A E = 106 cm -1 , 

k=^k\,: where — k\ = 1.05 x 10 -18 cm 3 mol _1 s _1 , 

(21) 

O2 ( 3 ^g)u=l +O2 ( 3 Zg) v= 0 2O2 ( 3 E g ) v= Q , 

k vT = k 3 = 400 . (22) 


method is distinctly better as compared to the analytical ap¬ 
proach for the present system. The reason being that the nu¬ 
merical solution requires solving simple algebraic functions 
as compared to solution of highly computationally intensive 
exponential and logarithmic functions in analytical approach. 
Another reason for preferring the former approach is the eas¬ 
ier prediction of variation in number densities at various loca¬ 
tions along and normal to the cell for a host of pumping possi¬ 
bilities. Figure 3 shows a typical case of the liquid oxygen cell, 
irradiated with a Gaussian pumping beam, discretized into 
a grid of 11 x 5. The numerical code employed has provisions 
for using wide ranges of gird arrays along with the possibility 
of employing non-uniform grids. This would allow packing of 
more number of grids in the region of larger gradients. Since 
the beam profile and the cell are symmetric with respect to the 
v-axis, the behavior of grids lying on the one side of the axis of 
symmetry will be identical with that on the other. Thus, only 
one half of the cell is considered for modeling to reduce the 
computational effort. 

In finite difference form, the rate equations (16) and (17) 
can be expressed as 

N u (t + A t) = rAt — k\N u (t)At — k 2 N u (t) 2 At — yN u (t)At , 

(24) 

N\(t + At) = rAt + ^ k\ N \x{t) At T yN u (t)At 

— k 3 N\{t)At. (25) 


The pump rate, r is governed by the following equation, 


r 



= noi 


A{hv) 


(23) 


where P\ is the incident pump power, 4 A’ is the focused beam 
area, 4 a ’ absorption coefficient (cm -1 ) corresponding to the 
pump wavelength , 4 n ’ is the number of ground state molecules 
pumped by each incident photon (n = 1 for 1064 nm, n = 2 for 
634 nm) to the singlet delta state, 4 hv’ is the incident single 
photon energy. The same relation is utilized for determining 
the rate of pumping to the first vibrational state as well, but is 
zero for the case of pumping wavelength being 634 nm and is 
a non-zero quantity for 1064 nm. 

Equations (16) and (17) are solved using finite differ¬ 
ence method in an explicit formulation. The finite difference 


For all cases in this study, the time step (At) for the solu¬ 
tion of the above equations is considered to be lower than the 
time associated with the fastest occurring physical process. 
For example, in case of continuous pumping the fastest oc¬ 
curring process is the decay of the upper lasing level, which 
is of the order of several tens of microseconds. Thus, in con¬ 
tinuous pumping case the assumed time step is in between 
1-10 |xs. At the same time, for the pulsed pumping case the 
fastest phenomenon is the pumping process itself and hence 
the time step duration is taken to be smaller than the pump 
laser pulse duration. The results presented in this paper have 
also been checked for independency of time steps. Further¬ 
more, in order to assess the cumulative response of the entire 
medium in the liquid cell, the representative densities of both 
levels at any given moment are obtained by spatially averag¬ 
ing the same over the generated grid. The relations employed 




Hr) = l 0 e 
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FIGURE 3 Grid geometry of optical 
pumping used for numerical modeling 
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for spatial averaging are given below, 


Mim if) — 


E n hj (0 


Mm(0 = 


l XJ 

E N kj (0 


i xj 


(26) 

(27) 


where 7’ and ‘ j ’ are the number of grid points in the V and 
V’ direction, respectively. Thus, the developed code not only 
predicts the response to the input pumping power at various 
locations along and across the cell but also predicts the aver¬ 
age response of the medium. It can be utilized for subsequent 
estimation of the range of expected gain and the correspond¬ 
ing extractable power output. 


4 Results and discussions 

For pure natural liquid oxygen, the upper laser 
level has a non-radiative lifetime (r u ) of about 44 |xs, while 
the lower level (r\) has relatively longer non-radiative life¬ 
time of about 2.5 ms [13]. This property of pure natural li¬ 
quid oxygen makes pumping very critical which is associ¬ 
ated with accumulation of population at the lower laser level. 
The rate of accumulation will determine the limit of popu¬ 
lation inversion duration and hence the small signal gain in 
this medium. This means, for a continuous population in¬ 
version and gain, that it is essential to either increase r u 
or decrease X\ or simultaneously do both by some means. 
Experiments have also proven that the lifetime of the sing¬ 
let delta oxygen is more in case of isotopic liquid oxygen 
medium as compared to that in pure natural liquid oxygen [14] 
medium. On the other hand there is experimental evidence 
that some impurity addition can play a role in the VT re¬ 
laxation rate of the lower level [ 8 ]. By considering these as¬ 
pects, three different oxygen-rich mediums in liquid phase 
viz: (i) natural liquid oxygen (ii) isotopic liquid oxygen and 
(iii) liquid air are considered in this study. We have numer¬ 
ically analyzed the evolution of population in both lasing 
levels and verified the possible population inversion dura¬ 
tion under different pumping conditions viz. wavelength and 
pump intensities. 

Two different wavelength sources at 634 nm and 1064 nm 
are considered for the analysis. Liquid oxygen has high ab¬ 
sorption coefficients for the above wavelengths, and hence it 
is not useful to have a gain length longer than 1 or 2 cm along 
the pump direction. Therefore, in the present study, a liquid 
cell of 1 cm length has been taken with the consideration that 
pumping beam of different power levels is focused to a fixed 
diameter. For the studies presented here, the beam energy is 
considered to be focused to 3 mm diameter and is passed only 
once through the length of the liquid cell. Simulations are 
carried out for various energies of the pump lasers to have 
various power densities and hence pump rates. The rate con¬ 
stants for various processes used in this simulation studies 
are mostly experimentally accepted ones, making the simu¬ 
lation results reliable. The time evolution of the population 
at the upper and lower levels obtained from (24) and (25) is 
averaged over the illuminated volume [(26) and (27)] and uti¬ 
lized to calculate the corresponding gain profile employing 


equation (14). The results obtained in each case are discussed 
below. 

4.1 Natural liquid oxygen 

In pure natural liquid oxygen, there are various re¬ 
ported values for singlet oxygen quenching with ground state 
molecules (k\) and its self-pooling rate (& 2 )- The reported 
values of self-pooling rate of singlet oxygen in liquid oxy¬ 
gen have wide conflicts. The earlier work by H. Klingshern 
et al. [14] had reported a value of 3 =b 1 x 10 “ 15 cm 3 /mol/s, 
while the recent experiments by K.R. Hobbs et al. [15] re¬ 
ported a value of 1 ±0.5 x 10 -17 cm 3 /mol/s, which was al¬ 
most a two-order difference. While the variation in quench¬ 
ing rates are minimal, we have considered the most accepted 
value of 2.2 x 10 4 s -1 , which is experimentally proved from 
the fluorescence decay signal [14]. The most accepted de¬ 
cay rate of the lower lasing level (£ 3 ), estimated theoretic¬ 
ally and observed experimentally, is 400 s -1 [7]. Therefore, 
the rate constants used for the simulation of natural liquid 
oxygen are k\ — 2.2 x 10 4 s -1 , & 2 = 1.0 x 10 -17 cm 3 /mol/s, 
k 3 = 400 s -1 and y = 0.25 s -1 . 

The simulation results for continuous wave (cw) pumping 
of natural liquid oxygen with 634 nm wavelength source for 
three different power levels of 1 W, 1 kW and 1 MWare pre¬ 
sented in Fig. 4a-c, respectively. The pumping of the same 
medium with 1064 nm source having 1 kW power is shown in 
Fig. 4d. The validation of the results is carried out by com¬ 
paring the singlet delta population in our estimations with 
the same obtained by the recently reported modeling results 
of K.R. Hobbs et al. [15], which they have also verified ex¬ 
perimentally. From Fig. 4b and d, it is clear that the gain in 
case of 634 nm source pumping is almost 35% more com¬ 
pared to that with 1064 nm source. It is clear that a positive 
gain is obtained in this medium for about 65 |xs and 57 |xs for 
634 nm and 1064 nm sources, respectively, is almost indepen¬ 
dent of the pump power. The singlet delta concentration and 
small signal gains are found to be almost linear with the pump 
power, which shows the possibility of near linear power scale 
up of the laser with the pump power. The deviations at very 
high pump intensities (~ MW/cm 2 ) are due to the dominance 
of pooling process at large singlet oxygen densities showing 
that efficient utilization of the medium needs optimized pump 
intensity. 

The results for pumping of natural liquid oxygen with 
10 ns pulsed laser at 1064 nm wavelength for two different 
energies of 1 J and 10 J are given in Fig. 5a and b, respec¬ 
tively. It is evident from the results that in case of pulsed 
pumping, obtained peak gain is linearly proportional to the 
pump energy for the reason that the decay processes are not 
able respond with in this time scale. The population inversion 
duration is found to be much larger than the pump pulse du¬ 
ration and hence natural liquid oxygen is a potential medium 
for pulsed lasers. The small signal gain obtained with 1 MW 
power cw pumping (0.02 cm -1 ) and 10 J, 10 ns pulsed pump¬ 
ing ( 0.01 cm -1 ) looks reasonable for obtaining lasing with 
better extraction efficiencies, provided the cavity losses can 
be minimized to less than 0.5%. These issues are addressed 
in forthcoming section on power extraction. It means that 
the medium can be utilized for lasing, if it can be pumped 
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Natural LO: 634nm, 1W Natural LO:634nm,1kW 





FIGURE 4 (a) The populations of the upper and lower lasing levels and the small signal gain for pumping with 634 nm source with 1 W power; (b) The 

populations of the upper and lower lasing levels and the small signal gain for pumping with 634 nm source with 1 kW power; (c) The populations of the upper 
and lower lasing levels and the small signal gain for pumping with 634 nm source with 1 MW power; (d) The populations of the upper and lower lasing levels 
and the small signal gain for pumping with 1064 nm source with 1 kW power 


Natural LO: 1064nm, 1J,10ns 


Natural LO: 1064nm, 10J,10ns 
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FIGURE 5 (a) The populations of the upper and lower lasing levels and the small signal gain for pumping with 1064 nm source with E = 1.0 J and 10 ns 

duration; (b) The populations of the upper and lower lasing levels and the small signal gain for pumping with 1064 nm source with E = 10.0 J and 10 ns 
duration 


with intensities of the order of few MW/cm 2 . But the only 
concern with natural liquid oxygen is the time duration of 
positive gain, which limits its use as a possible cw lasing 
medium. 

4.2 Isotopic liquid oxygen 

The lifetime, r u , of C> 2 ( 1 A g ) i;= o is experimentally 
found to be improved in isotopic liquid oxygen ( 18 02 ) from 


44 pis to about 10 ms [14]. The quenching of 02( l A g ) v =o 
molecules with the ground state molecules O 2 ( 3 ^) i; =o has 
a drastic isotopic effect, making k\ considerably smaller in 
isotopic liquid oxygen 18 02 than for natural liquid oxygen 
16 C> 2 . The reasons for this are detailed in [13]. The above work 
also reports that the V-T relaxation rate for the first vibra¬ 
tional level of ground molecules in isotopic liquid oxygen is 
still 400 s -1 , which is the same as that for natural liquid oxy¬ 
gen. At the same time, H. Klingshim et.al have reported that 
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Isotopic LO: 634nm, 1 kW 


Isotopic LO: 634nm, 1 MW 
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Isotopic LO: 1064nm, 1 kW 
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FIGURE 6 (a) The populations of the upper and lower lasing levels and the small signal gain for cw pumping with 634 nm source with 1 kW power; (b) the 

populations of the upper and lower lasing levels and the small signal gain for cw pumping with 634 nm source with 1 MW power; (c) The populations of the 
upper and lower lasing levels and the small signal gain for cw pumping with 1064 nm source with 1 kW power; (d) The populations of the upper and lower 
lasing levels and the small signal gain for pulsed pumping with 1064 nm source 10 J, 10 ns pulse 


the pooling rate of C> 2 ( 1 A g ) i;= o molecules in isotopic liquid 
oxygen is an order of magnitude higher than in natural liquid 
oxygen [14]. 

Based on the above discussions, the rate constants con¬ 
sidered for the analysis on isotopic liquid oxygen are the 
following: k\ = 100s _1 , = l.Ox 10 -16 cm 3 /mol/s, k 3 = 

400 s -1 and y = 0.25 s -1 . The simulation results for con¬ 
tinuous wave pumping with 634 nm wavelength source for 
two different power levels of 1 kW and 1 MW are shown in 
Fig. 6 a and b, while the results for pumping with 1064 nm 
source for 1 kW cw power and 10 J, 10 ns pulse are given 
in Fig. 6 c and d respectively. The results indicate that pump¬ 
ing with 634 nm source yields a continuous population inver¬ 
sion for very large duration. At the same time while pump¬ 
ing with 1064 nm source, the population inversion duration 
is reduced to only 10 ms, which is still larger compared to 
that of natural liquid oxygen. However, as in the case of 
natural liquid oxygen, the population inversion duration is 
found to be independent of the pumping power. Further, 
for low pump power (for example the case of 1 kW), it is 
found that the small signal gain in isotopic liquid oxygen 
is large as compared to that in natural liquid oxygen. How¬ 
ever, the small signal gain is almost the same for the case 
of high pump power (~ 1 MW). The reason for this result 
is that at very high pumping intensities, pooling becomes 
the dominant loss mechanism compared to that of quench¬ 


ing, which will balance the population in both levels much 
faster. 

4.3 Liquid air 

In case of natural liquid oxygen and isotopic liquid 
oxygen discussed earlier, the lower level lifetime remains the 
same, while the upper level lifetime is found to improve for 
the latter. In these studies, the lifetime of the lower lasing 
level mainly depends on the vibrational to translational con¬ 
version (V-T) and is about 2.5 ms. But there is clear evidence 
from the experiments that the V-T relaxation rate constants 
are highly dependent on the impurities (like Ar, Kr, Xe, CO 2 , 
H 2 O etc.) present in the liquid oxygen, and experimentally it 
was observed to be ~ 50 |xs by Renner et al. [ 8 ]. This dras¬ 
tic reduction of the lower level lifetime by the impurities may 
be explored and gives hope for a reasonably long continu¬ 
ous population inversion condition in impure liquid oxygen. 
The influence of the impurities on the lifetime of the ( l A g ) v= Q 
molecules has not been reported by any one to the best of our 
knowledge. Furthermore, it is also interesting to note that the 
quenching rate of the upper level is found to reduce almost 
linearly with the percentage of oxygen present in the mixture 
along with buffer medium such as O 2 /N 2 or 02 /Ar etc., which 
improves the lifetime r u [2] as well. However, it should also be 
remembered that the absorption coefficient for the pump beam 
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Liquid Air: 634nm, 1 MW Liquid Air, 1064nm,10J,10ns 



FIGURE 7 (a) The populations of the upper and lower lasing levels and the small signal gain for cw pumping of liquid air medium with 634 nm source hav¬ 

ing 1 MW power; (b) The populations of the upper and lower lasing levels and the small signal gain for pulsed pumping of liquid air medium with 1064 nm 
having 10 J, 10 ns 


will be reduced with oxygen density reduction [5], which may 
require multiple passage of the beam for efficient pumping. 

In order to explore this idea, liquid air, where there is 
availability of almost 21%-22% of liquid oxygen and rest 
liquid nitrogen has been considered for the analysis. Here 
we have assumed that the liquid air has sufficient impuri¬ 
ties to reduce the lifetime of the lower vibrational levels to 
about 50 |xs. Since nitrogen is a buffer medium, it typic¬ 
ally does not affect the lifetime of electronic energy levels. 
The reduction in the quenching and pooling rates is propor¬ 
tional to the oxygen density in liquid air, which is less as 
compared to that in pure liquid oxygen. Therefore, for the 
analysis of liquid air, we have considered rate constants as 
thefollwoing: k\ = 4.4 x 10 3 s _1 ,fc = 2x 10 -18 cm 3 /mol/s, 
&3 = 20000 s _1 and y = 0.25 s -1 . The studies are carried out 
for both cw and pulsed pumping for various intensities at 
634 nm and 1064 nm. However, the simulation results for the 
case of continuous wave (cw) pumping with 1 MW source 
at 634 nm and the same for pulsed pumping with 10 J, 10 ns 
pulse at 1064 nm are only presented here in Fig. 7a and b re¬ 
spectively. 

In case of pumping with a continuous source, it is found 
that a stable gain and continuous population inversion con¬ 
dition occurs while pumping with either of the wavelength 
sources unlike in other two discussed mediums. The gain ob¬ 
tained in liquid air at lower pumping power (say 1 kW or less) 
is relatively less compared to both natural and isotopic liquid 
oxygen due to the reduction in the excitation rate because of 
the low oxygen density available. But for very high pump 
power (~ 1 MW), the gain obtained with liquid air medium 
is much higher compared to both natural as well as isotopic 
liquid oxygen. This clearly indicates the role of self pooling 
at high density production of upper lasing level molecules 
in mediums characterized with high oxygen density (natural, 
isotopic oxygen). In case of pulsed pumping with nanosec¬ 
ond pulses, the gain is almost linear with the pump energy and 
nearly the same for all the three mediums. 


5 Output power estimation 

The extractable output power from the lasing 
medium is proportional to the available power in the medium 


and is expressed as: 

Fext — Pav^lex t > ( 28 ) 


where ^ ex t is the extraction efficiency expressed by the prod¬ 
uct of the medium extraction efficiency (rj ex tm) and the res¬ 
onator extraction efficiency (i/ ex tr)- The medium and resonator 
extraction efficiency can be estimated using the following set 
of Rigrod relations [16] as 


^extm — 


t 

7?e X tr — — 

t -\-a 


^ , 

YoJ 

_ y th (2L)-a 

Kth(2L) 


(29) 

(30) 


where the threshold gain (]/ t h) essential for lasing is dependent 
on the gain length L, resonator coupler mirror transmission (t) 
and the absorption loss ( a ) of the resonator mirrors, which can 
be expressed as 


t T ci 



(31) 


For mediums having very low stimulated emission cross- 
section, very large pumping power is essential to obtain the 
threshold gain to overcome the resonator losses. In such medi¬ 
ums, typically, the available gain will be near to the threshold 
gain and hence the medium extraction efficiency will be very 
low. However, if the medium has high saturation intensity, it 
can have high amplification through proper selection of the 
resonator mirrors and can yield high output power. The calcu¬ 
lation shows that the extraction efficiency for this medium is 
highly dependent on the small signal gain. Hence, for a high 
gain case of 0.06 cm -1 , an extraction efficiency of about 60% 
can easily be obtained provided the resonator loss is limited 
to less than 0.5%. The extractable energy in case of pulsed 
pumping and possible output power in case of continuous 
pumping are estimated employing the relations [10], 


£out = E^ATil-a) 
Pout — hatAT(l — a) 


-2ypL 
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FIGURE 8 (a) The output power with coupler transmission for natural liquid oxygen and liquid air under different gain conditions and continuous pumping; 

(b) Change in output power with coupler transmission for isotopic liquid oxygen under different gain conditions during continuous pumping 


where E sat , A are the saturation energy density, area of cross- 
section of the gain medium respectively and a is the absorp¬ 
tion loss in the medium. The saturation intensity for natural 
liquid oxygen and liquid air are considered to be the same 
due to their identical upper level lifetime. But the saturation 
intensity for the isotopic liquid oxygen will be much less be¬ 
cause of its relatively large lifetime and is estimated to be 
~ 230kW/cm 2 . 

In case of continuous pumping, the estimated output 
power for the gain values predicted by modeling ranging from 
0.005 cm -1 to 0.06 cm -1 corresponding to natural liquid oxy¬ 
gen and liquid air are plotted in Fig. 8a and the same for 
isotopic liquid oxygen is shown in Fig. 8b. For the observed 
gain of 0.06 cm -1 for 1 MW pumping in liquid air, the analy¬ 
sis predicts a photon conversion efficiency of more than 25%. 
At the same time, for pure natural liquid oxygen it is only less 
than 5% due to the smaller gain obtained in the medium for the 
same pump power. Further, the photon conversion efficiency 
is found to be too low in isotopic liquid oxygen, i.e., only 
0.02% compared to 4.5% in liquid air for the same gain value 
of 0.015 cm -1 for the reason that it has very low saturation 
intensity. Typically, the small signal gain is found to be a func¬ 
tion of the pump wavelength. Thus, the gain obtained is almost 


Energy vs.Coupler transmission for Nat. Lo,Isotopic LO, Liquid Air 



FIGURE 9 Change in output energy with coupler transmission for natu¬ 
ral, isotopic liquid oxygen and liquid air under different gain conditions and 
pulsed pumping 


35% more in case of 634 nm pumping as compared to that 
of 1064 nm pumping, the corresponding photon conversion 
efficiency is also found to be proportionally higher. 

In case of pulsed pumping, the energy output is almost 
the same for all mediums and dependent only on the gain 
since E sat is dependent only on the stimulated emission cross- 
section. The predicted energy output variation with a coupler 
transmission for different gain values is shown in Fig. 9. The 
results indicate that an energy conversion efficiency of about 
16% is possible with these mediums. 

6 Realization criticalities 

Major practical challenges in the realization of li¬ 
quid oxygen laser have been addressed along with certain 
solutions in [17]. Apart from these, one such criticality antici¬ 
pated is the thermal distortions of the beams due to the heating 
of the lasing medium. This occurs due to the chemical kinetics 
of the molecules in the medium during the excitation pro¬ 
cess. The V-T relaxation process (22) is the one, which results 
in slow heating of the solution with energy release of about 
0.19 eV in each transition. We have estimated the temperature 
rise by incorporating this energy release with the correspond¬ 
ing decay rate and found that nearly l%-3% of the pump 
energy is dumped resulting in heating of the medium for pump 
power of 1 MW to 1 kW. At high pump power the fraction 
is less because the pooling process dominates the other de¬ 
cay processes viz, V-V and V-T relaxations. In the case of 
pulsed pumping with low repetition rate, it is found that the li¬ 
quid medium will not undergo any physical change due to the 
small amount of energy dumped during the pump duration. 
In case of very high power (1 MW) cw pumping, it is found 
that the heat dumping can affect the physical stability since 
the heat dump rate is more than the latent heat of vaporization 
(2934 BTUlb _1 mol _1 ) of the liquid. In such cases, it will be 
essential to use either super-cooled liquid with cooling provi¬ 
sions to maintain its temperature or to use replenishing lasing 
medium to avoid the temperature rise of the medium. 

The temperature rise of the medium can also cause a re¬ 
fractive index profile within the medium and may induce 
thermal blooming effects and distortion of the beam. Experi¬ 
mental studies on such thermal blooming effects and beam 
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distortions in liquid nitrogen are discussed in [18]. In case of 
pulsed pumping with less repetitive pulses the natural con¬ 
vection will take care of the thermal effects hence reducing 
beam distortions. But in case of steady state cw pumping or 
high repetitive pulsed pumping, these distortions will be tran¬ 
sient for initial pump duration and will reach an equilibrium 
condition for the remaining duration. In the present studies 
we have considered the cases where around 80% (wavelength 
dependent) of the beam is absorbed within the medium and 
estimated the average populations of the various levels and 
hence will not affect the computed results. 

7 Conclusion 

The present studies indicate that oxygen rich cryo- 
solutions can be an effective lasing medium with excellent 
photon conversion efficiency. The results indicate that rea¬ 
sonable small signal gain for lasing action is achievable only 
when pump intensity is close to MW/cm 2 . The simulations in¬ 
dicate that in all the three mediums considered here, pumping 
with 634 nm source gives higher small signal gain in both con¬ 
tinuous and pulsed mode as compared to that with 1064 nm 
source. The analysis shows that liquid air medium is highly 
potential for continuous wave lasers when pumping with ei¬ 
ther 1064 nm or 634 nm source while isotopic liquid oxygen 
is suitable only when pumping with 634 nm source. It is also 
found that pure natural liquid oxygen is a suitable medium 
only for pulsed lasers and continuous gain is available only 
for few milliseconds. Liquid air shows reasonable small sig¬ 
nal gain, and for a case of MW level pumping, it is about 
0.06 cm -1 with a corresponding photon conversion efficiency 
of about 25%. Although, isotopic liquid oxygen is suitable 
for continuous wave lasing, it provides a low photon conver¬ 


sion efficiency because of the low saturation intensity of the 
medium. Also, for very large pump intensity, pooling loss is 
the dominant process governing the small signal gain. Hence, 
liquid air is the most potential medium, since it is associated 
with lower pooling losses due to the presence of buffer nitro¬ 
gen. We are in progress of setting up the experimental setup 
for the studies with these mediums, which will surely shed 
more light on this possible lasing medium. 
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